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Abstract 

We look for physically realistic initial data in numerical relativity which are in agreement with 

(N 

' post-Newtonian approximations. We propose a particular solution of the time-symmetric con- 

o : 

■ straint equation, appropriate to two momentarily static black holes, in the form of a conformal 

o 
m 
o 



decomposition of the spatial metric. This solution is isometric to the post-Newtonian metric up 
to the 2PN order. It represents a non-linear deformation of the solution of Brill and Lindquist, 



^ • i.e. an asymptotically flat region is connected to two asymptotically flat (in a certain weak sense) 

5h . sheets, that are the images of the two singularities through appropriate inversion transformations, 

bp. 

The total ADM mass M as well as the individual masses m-i and m2 (when they exist) are com- 
^ ' puted by surface integrals performed at infinity. Using second order perturbation theory on the 

Brill-Lindquist background, we prove that the binary's interacting mass-energy M — mi — m2 is 
well-defined at the 2PN order and in agreement with the known post-Newtonian result. 
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I. MOTIVATION AND RELATION TO OTHER WORKS 



The numerical computation of the colhsion of two black holes is of paramount importance 
for the observation of gravitational waves by the network of laser-interferometric detectors. 
When investigating this problem the ten Einstein field equations are separated into: (i) four 
constraint equations that are to be satisfied by some initial data given on an initial three- 
dimensional Cauchy hypersurface; (ii) six hyperbolic-like equations describing the dynamical 
evolution of the gravitational field on neighbouring hypersurfaces. The Bianchi identities 
guarantee that the constraint equations are satisfied on neighbouring hypersurfaces if they 
are on the initial hypersurface. There are infinitely many ways that the initial data can be 
chosen to represent the starting state of the evolution of black holes. It is widely admitted 
that the problem of choosing physically realistic initial conditions for the collision of two 



black holes has not yet been solved. There has been a lot of concern in the litterature 
for knowing what would really motivate physically a particular choice of initial data. 

Let us consider the problem of time- symmetric initial data, which are physically appro- 
priate to two momentarily static black holes, i.e. with zero initial velocities. The dynamical 
evolution of time-symmetric data describes the subsequent head-on collision of the two black 
holes. In this situation the constraint equations reduce to the Hamiltonian or scalar con- 
traint equation i? = (in vacuum — the case appropriate to black holes), with R being the 
three-dimensional scalar curvature. Considering as usual 0,0,0] a conformal decomposition 
of the spatial metric (spatial indices i,j,--- = 1, 2, 3), 



where 'jij is the physical metric and denotes the conformal (unphysical) metric, we 
obtain the Lichnerowicz 0] equation, which is an elliptic-type equation to be satisfied by 
the conformal factor \E'. In the time-symmetric case that equation becomes 



The scalar curvature R and Laplacian A are the ones associated with the conformal metric 
7jj. The interest of the conformal decomposition is that by solving Eq. ()1.2j) we generate a 
physical solution of the constraint equation starting from any choice for the conformal 
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metric 7^-, (a priori). Therefore the problem of initial conditions resides in the choice of a 
physically well-motivated conformal metric jij. 

The simplest choice of initial conditions (one motivated by simplicity rather than by 
physics) is the one for which jij equals the fiat metric 6ij. In that case Eq. reduces to 
the flat-space Laplace equation. Some exact solutions, appropriate to (momentarily static) 
)lack holes, have been obtained by Misner 

7|. The solution of Brill and Lindquist [for which the conformal factor takes the form 
of Eq. ()3.2b|) below] is particularly interesting: it describes the "geometrostatics" of two 
black holes, consisting of three asymptotically flat regions connected by two Einstein-Rosen 
bridges (actually, the solution is known for black holes). The one region containing the 
two throats is supposed to represent our universe, while the two sheets expanding behind 
it are associated with the two black holes. The beauty of the Brill-Lindquist solution is 
that not only the total ADM mass-energy M of the binary, but also the two individual 
masses mi and m2 of the black holes, are computed "at infinity" . One can use for instance 
standard surface integrals extending on topological two-spheres at infinity. The binary's 
geometrostatic energy (i.e. the gravitational interacting or binding energy, in the center-of- 
mass frame) is therefore computed unambiguously as ^ = M — mi — m2. 

Yet the solution of Brill and Lindquist, despite its undeniable interest, is not "physically 
realistic" in the sense that it differs, in the limit c —>■ +00, from the three- metric found 
for a post-Newtonian solution. Indeed the post-Newtonian metric generated by two point- 
particles is known to deviate from conformal flatness at the 2PN order (see e.g. Ref. |l2]) ^■ 
In consequence both the Brill-Lindquist metric and the associated binding energy E disagree 
with the post-Newtonian results from the 2PN order. 

A compelling motivation for constructing physically realistic initial data is the agreement 
with the post-Newtonian approximation when c — +00. Recall that the post-Newtonian 
theory provides some explicit expressions for the metric, equations of motion and energy of 
binary systems of point-particles. The post-Newtonian metric is valid in the binary's "near- 
zone" (size of near- zone is much less than a gravitational wavelength), but has been proved to 
come from the re-expansion when c — > +00 of a "global" (post-Minkowskian-type) solution. 



^ In this paper c and G denote the speed of hght and the gravitational constant. As usual the nPN order 
means the terms of order when c — > +00. 
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defined everywhere in space-time including the wave zone Furthermore, in the post- 
Newtonian approach the modelling of black holes by point-like particles — i.e. technically 
by Dirac delta-functions in the stress-energy tensor — is rather well justified. We shall 
provide below some further evidence that the "post-Newtonian masses" of point-particles 
are indeed identical to the black- hole masses^. 

It has been suggested [3| (see also 17 1) that, in order of taking into account the post- 



Newtonian physics into the initial data, one should adopt for the conformal metric jij 
directly a post-Newtonian solution. In such a proposal one expects that by correcting the 
post-Newtonian solution by means of a conformal factor (computed numerically), one 
will somehow be able to "compensate" in the physical metric 'jij (i.e., in fact, to cancel out) 
the systematic higher-order post-Newtonian error terms that are neglected in jij. 

There has been other proposals for realistic initial data, in particular built on the re- 



laxation of the assumption of conformal flatness One of these is to use for 'jij a linear 
combination of two (boosted versions of) Kerr-Schild metrics 3,^^- It is not known if the 
physical metric which is generated in this way from the numerically-computed conformal 
factor, is consistent with post-Newtonian (e.g. 2PN) calculations. 

On the other hand, we should remark that to which extent the hypothesis of conformal 
flatness introduces some unphysical spurious (and numerically important) effects remains 
an issue. Indeed, relaxing this hypothesis may not always be a panacea. A quite different 
idea for settling the initial conditions of two black holes is to solve numerically a subset 
of the Einstein field equations (the four constraint equations plus one evolution equation), 
under the two premices of conformal flatness and the existence of a helical Killing symmetry 



20 



21[. Despite the conformal flatness of the spatial metric, the latter calculation gives 
results in good agreement with post-Newtonian predictions 0,0, Is^. Nevertheless, we 
think that it is important, at some stage, to get rid of the hypothesis of conformal flatness. 
We emphasize furthermore that the agreement between numerical relativity 0, ^ and 



post-Newtonian theory 
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holds up to the very relativistic regime of the innermost 



circular orbit (ICO), where the orbital velocities are of the order of 50% of the speed of 



light ^. In Refs. 



22| it is suggested that the result for the ICO of two black holes with 



^ See also, in a similar context, the matching of a IPN solution for the orbital motion of point-particles to 

two perturbed Schwarzschild black holes [l5j . 
^ The ICO is defined by the minimum of the binary's energy function for circular orbits. It represents 
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comparable masses, at the 3PN approximation^, is likely to be close to the "exact" solution, 
within 1% of fractional accuracy or better. This constitutes a motivation for advocating 
that the black-hole initial conditions, which are to be set (presumably) around the location 
of the ICO, should be in agreement with post-Newtonian theory. 

In the present paper we propose an alternative way for incorporating the post-Newtonian 
information into the initial data of black- hole binaries (in the time-symmetric case). We 
find, in Section |nj a simple expression for the conformal metric jij, which is such that the 
corresponding physical metric 7,^ is isometric (i.e. differs by a coordinate transformation) to 
the standard post-Newtonian spatial metric at the 2PN order in the limit c +00. At the 
same time, the solution is defined globally in space, with a global structure similar to the 
one of Brill and Lindquist. Our solution is not "exact", but exists as a certain non-linear 
perturbation, investigated in Section lTTTl of the Brill-Lindquist solution, playing here the role 
of a "background" metric. Most importantly, in Section IIVI we investigate the asymptotic 
structure of the solution, and compute "geometrically" the masses M and mi, m2, i.e. by 
surface integrals performed in their respective domains at infinity. In Section |3 the binary's 
interacting energy, deduced from the previous masses, is proved to be in agreement with 
the known post-Newtonian energy up to the 2PN order. (We shall find, however, that the 
definition we adopt for the two individual masses mi and 1712 makes sense only up to the 
2PN order.) 

II. CONFORMAL DECOMPOSITION OF THE SPATIAL METRIC 
A. Definition of the conformal metric 

The conformal metric we propose, appropriate to two black holes of masses mi and 1712 
(to agree later with the "geometrical" masses), located at the singular points yi and y2, and 

a useful reference point for the numerical ji^, 13 ^^"^ 3PN 22] calculations because for both of them 
the ICO is well-defined and can be meaningfully compared. However, the radiation reaction terms are 
neglected in its definition, so the ICO probably does not have a rigorous physical meaning in a context of 
exact radiative solutions. 

^ We mean the standard Taylor-expanded form of the approximation — without using any post-Newtonian 
resummation techniques. 
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momentarily at rest (vi = V2 = 0), is 



= — 7^:::^^ ■ ^^.i) 



The function g introduced here represents an elementary "kernel" playing an important role 
in post-Newtonian calculations . It depends on the "field" point x on 

the one hand, and on the pair of "source" points yi, y2 on the other hand; it is defined by 



fl'(x;yi,y2) = ln(ri + r2 + ri2) , (2.2) 

where ri = |x — yi| and r2 = |x — y2| are the distances to the black holes, and ri2 = |yi — y2| 
is their separation. The function g satisfies, in the sense of distributions (i.e. for all values 
of X, including the singular points yi and y2), 

Ag = ^, (2.3) 

rir2 

where A denotes the usual fiat-space Laplacian with respect to the field point x. In Eq. 
()2.1|) the derivatives are taken with respect to the two source points^. The carets around 
the indices refer to the symmetric and trace-free (STF) projection: T<jj> = ^ (Tij + Tji) — 
l^ijTkk] so the trace of the metric ()2.1|) is normalized to be 7,^ = 3. It is of course quite 
natural (and useful in practice) to impose that the deviation of the conformal metric from 
fiat space be trace-free. 

Our proposal is to generate, by means of numerical techniques (i.e. with the help of 
elliptic solvers), a conformal factor \I' solving the constraint equation ()1.2j) corresponding to 
the particular choice of conformal metric ()2.1|) . The metric = \l/^7jj we obtain in this 
way will incorporate the post-Newtonian (2PN) physics of the initial state of the head-on 
collision of two black holes ^. The latter assertion will now be proved, for the rest of the 
paper, with the help of analytic perturbation methods. 
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The following explicit formula holds: 



^ d^g ^ »l2»i2 - ^ (nia -ni)(nj2 +n^) 

'^■^ dyidyi ruin + r2 + ri2) (ri + r2 + ^12)2 

Here, n\ = (x' — y{)/ri, n\ = (x* iPiil'^'i '^12 = (j/i ^ y2)/''i2- See Ref. I^l for further discussion 
and formulas, CliD-cerning the function g. 
6 See Refs 



24, l25l| for numerical calculations of the head-on collision of black holes; see also Ref. I2a for 



a post-Newtonian calculation. 



6 



B. Relation with the post-Newtonian metric 



The first result shows that the "near- zone" behaviour of the solution (i.e. r/c — > 0) is 
physically sound. 

Theorem 1 The post- Newtonian expansion (when c — > +00) of the solution = ^'^^j of 
the constraint equation 1^1.^ . where is defined by Eq. {\2.1\ . differs from the standard 
post-Newtonian spatial metric, calculated by standard post- Newtonian methods, by a mere 
change of coordinates at the 2PN order, i.e. 



7.,=gri +^^i,+^,i. + 0[-] . (2.4) 



lvi=V2 = \ ^ 



1 



To be precise, by gff ^ we mean the spatial metric in harmonic coordinates, when trun- 



3«J 



cated at the 2PN order, that is given by Eq. (7.2c) in Ref. As indicated in Eq. ()2.4j] . 
we must set the particles' velocities vi and V2 to zero in the post-Newtonian metric in order 
to conform with the assumption of time-symmetry. The remainder O stands for the 
neglected 3PN and higher-order terms. 

The proof of Theorem Q is easily achieved on the basis of a post-Newtonian iteration of 
Eq. (|1.2|) . At the 2PN order, this equation becomes 

A* = -^^9, (<*>) + O (i;) . (2.5) 

where we denote iOj = ^ ^ f , . The most general solution reads as 

* = V'-^D(f + ^Vof^l . (2.6) 



2c^ \3 2ri2 J Vc^ 
where ip represents a solution of the homogeneous equation (i.e. Aip = 0) ^. By comparing 



with the post-Newtonian metric (see Eq. (7.2c) in Ref. lj| in which vi = V2 = 0), we 



^ We employ the notational shorthand D = gjigyi (so that Dg = ,3^). Notice the useful relations 

X. 1 1 1 

D5 



2rir2 2riri2 2r2ri2 
1 1 1 



2rir2 2riri2 2r2ri2 



7 



readily find that the latter solution is uniquely specified, at the 2PN order, as being 



Grrii ( Gmo \ Gnio ( Grrii \ / 1 \ , , 

Also uniquely determined is the expression of the vector in Eq. ()2.4|1 . which represents 
an infinitesimal gauge transformation. At the 2PN order it takes the expression 

where rii = (x — yi)/ri and n2 = (x — y2)/T2- The system of spatial coordinates employed 
in the present paper, i.e. which corresponds to a conformal metric of the particular form 
displayed by Eq. ()2.H) . is thereby determined, and related up to the 2PN order to the 
harmonic coordinate system used in Ref. by the coordinate change = ^' (with 

Summarizing this section, we have found a conformal decomposition of the metric 
that is in agreement with the post-Newtonian metric up to the 2PN order. Notice, however, 
that the conformal metric we propose is not unique, because we can always add higher-order 
terms (e.g. 3PN) to Eq. ()2.H1 . A contrario, this means that with the metric ()2.H1 one should 
not expect Theorem to hold at the 3PN order and beyond. At the 3PN order for instance 
the conformal metric will be more complicated than the simple expression ()2.1|) . 

III. PERTURBATION OF THE BRILL-LINDQUIST SOLUTION 

Having checked the near-zone structure of the metric, let us investigate some of its global 
properties, when it is viewed as a solution of the constraint equation p.2j) that is a priori 
valid everywhere on a space-like hypersurface. For this purpose, we impose that 'jij is a non- 
linear perturbation of the exact solution of Brill and Lindquist 0] — which is conformally 
fiat {•y^^ = 6ij). This will mean that the topology of our solution is identical to the topology 
of the Brill-Lindquist solution, i.e. be "three-sheeted", in contrast with the two-sheeted 
topology of the Misner-Lindquist solution ^, . 

A. Hierarchy of perturbation equations 

Let us write Eq. ()2.1|) in a more transparent form, = 6ij + /Ssij, where 
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/? = , (3.1a) 

Si~ = : ^ (such that Su = 0) . (3.1b) 

It is helpful to view the non-linear perturbation we want to consider, as being generated by 
the "seed" or "generating" function Sij, and to interpret the parameter (3 as the magnitude 
of that perturbation. In the following, j3 will play the role of a "book-keeping" parameter al- 
lowing us to label the successive non-linear perturbation orders. Besides the non-conformally 
fiat piece of the metric brought about by Sij, it is evident that we must also introduce a 
perturbation in the conformal factor. We pose 



^ = iJ + a , (3.2a) 

^ = 1 + ^ + ^, (3.2b) 
ri r2 

where ip denotes the Brill-Lindquist conformal factor parametrized by two constants ai 
and a2 (in the case of two particles). The quantity a denotes a certain perturbation of ip. 

With full generality — within the present perturbative framework — , we look for the 
expression of cr in the form of an infinite power series in (3, solving the equation ()1.2|) : 

+ 0O 



a = J]/5V(„). (3.3) 



n=l 



We insist that in our approach the Brill-Lindquist solution plays the role of the background 
metric (it depends solely on ai and 02), while the function Sij yields a non-linear deformation 
of this background, perturbatively ordered by the book-keeping parameter (3. For a given 
choice of Sij, we expect that the resulting solution is unique (at least in a sense of formal 
power series in /?). Because f3 involves two mass factors mi and m2, so it is proportional to 
G^, the perturbation series we look for is like a "double" post-Minkowskian expansion, going 
"twice as fast" as the usual post-Minkowskian expansion when G — 0. Considering this 
series on the point of view of a post-Newtonian re-expansion, we see that each non-linear 
order brings in a new factor l/c^, so that our perturbation series can be said to go by — 
quite efficient indeed — steps of 2PN orders. 
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Let us compare our definitions ()3.2|) - ()3.3p with the result ()2.6|) - ()2.7|) provided by the 
agreement with the 2PN metric. We find that the constants ai,a2 are determined with 
relative IPN accuracy, 



GlTLi 



and that with this accuracy they agree with the prediction of the Brill- Lindquist solution^. 
In addition, we find that the perturbation a in the conformal factor is given by the second 
term in the right side of Eq. ()2.6|) . which comes in only at the 2PN order — it is purely 
linear in j3. Hence, 

with all higher-order cr(„)'s being negligible with this approximation. Actually, Eqs. ()3.4|) - 
fl3.5|) correspond to a particular "sharing" of the terms between ai, 02 on the one hand, and 
cr(i) on the other hand, since we can always add to cr(i) some "homogeneous" terms ~ 1/ri 
and ~ l/r2 without changing the equation for (T(i). It is obvious that such a sharing is 
physically irrelevant. However we shall forbid a different sharing of terms by adopting below 
the prescription (which represents simply a convenient choice) that cr(i) and all subsequent 
iterations solve the equation for the conformal factor in the sense of distributions. Some 
results more complete than ()3.4j) - ()3.5|) will be obtained below. 

By inserting the perturbation ansatz ()3.3j) into the constraint equation ()1.2|1 . and by 
identifying each of the coefficients of the successive powers of (3 in both sides of the equation, 
we obtain a hierarchy (indexed by n G N) of Poisson-type equations: 

Aa(„) = [cr(i), ■ ■ ■ , a(„_i)] , (3.6) 
where we recall that A = didi. The nth-order source term S(„) depends on the solutions 



We recall that and af^ in the case of the Brill-Lindquist solution are related to the masses by the 
exact relations 



BL 



1 + ^- 

ri2 



„BL I „BL 
"1 +"2 - ^ 



— and 1^2 

GM 



10 



of the preceding iterations (T(i),-- - ,(J(„_i), on the "generating" function Sij and on the 
"background" conformal factor ip. Notice that ip satisfies the equation Atp = — 47r(ai5i + 
02^2), where e.g. 61 = 5(x — yi) denotes the Dirac function at the point yi. Obviously we 
may include the equation for the background conformal factor into our hierarchy of equations 
by posing ^ = Y.n=o f^'^'^in)^ with (j(o) = ip and E(o) = -An{ai6i + a2S2). 



B. Analytic closed form of the linearized solution 

To the linearized order the perturbation equation reads explicitly 

Acr(i) = ^i^dijSij + djSijdiij + Sijdijip . (3.7) 

o 

Interestingly, this equation turns out to be solvable in analytic closed form. We find that 
the unique solution of Eq. ()3.7p . that is valid in the sense of distributions and tends to zero 
at spatial infinity (i.e. when r = |x| ^ +00), is 



^(1) 



16 

+ ai 



-D 



ri + r2 



12 



2 32 32 



_9_ 

ri2 



32 





In 


'n - 






■ri2- 






A- 







9 

+ — D 
32 



+ 



Inri 

ri2 
1 



+ 



^9 
12ri 



2Ariri2 32r2rf2 



a2 



{l„2} 



(3.8) 



Besides the already met shorthand D 



source points by Ai 



92 



dy\dy\ 



and Ao 



, we denote the Laplacians with respect to the 
QyiQyi ■ In Eq. fj3.8|) . the two terms proportional to 

ution involves the 



dy\dyi^ 
92 
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12| | for solving 



ai and 02 are deduced from each other by label exchange 1^2. The so 
special functions Hi and Ki (and 1 2) which were introduced in Refs. 
some elementary equations (in the sense of distributions) in the problems of equations of 
motion and wave generation at the 2PN order. We have, 



AiJi = 2igjdij 
AfsTi = 2 0^ 



Inri 



r2 



(3.9a) 
(3.9b) 
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in which jfif, = ^ ^ f , . The functions Hi and Ki admit the closed- form expressions given 
by Eqs. (3.48)-(3.51) in Ref. or equivalent^ by Eqs. (6.3)-(6.5) in Ref. tl^J; for the 
present purpose we adopt the formulas^ 



Hi = Ai 



2ri2 



r2 



+ 



2rlrl2 2r'(ri2 2rir^2 



Ki 



D ( — In 

T2 



ri 



1 



D 



1 



lnri2 

ri 



2 V ru 



r2 



ri2-l / 2r^r2 2r2ri2 2rfr 



12 



(3.10a) 
(3.10b) 



Anyway, we find that the exact solution of the linearized perturbation equation is 



(1) 



(- + ^^^^^ 
16 \3 2ri2 



_9_ 

ri2 



-1a, (a 

32 -\ri2 



15r2 



2 V ^1 

1 



QAr'lr'l2 4rjri2 64rjr2 24ri 



--d( 


^lnri\ _^ 




32 \ 




12ri 


7 


3 


} 


24rir^2 


64r2r^2 





+ a2|l 



(3.11) 



Let us quote also, for completeness, the fully explicit form obtained by expanding all the 
derivatives in this result: 



1 



^(1) 



ri + r2 



96rir2 64rf2 192ri2 v^i 



1 1 

^2 



1 



— + — 

64r5*2 V'"2 ^1 



+ "1 < --r^ - 
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+ 



ri 



+ 



4rf 64rf2 24rirf2 192r^ri2 192rfr2 32r2rf2 64r2r^2 



1 



+ 



^2 



ISro 



+ 



^^2 



+ 



15r2 



+ 



24rir2ri2 4rir^2 64r^r^2 4rjri2 64r^r^2 4rjr^2 



4rjrj2 



® We notice here (since it was not noticed in Refs. [loi H^ ) that Ki can also be given an interesting form 
in connection with a simple elementary kernel fci, viz 

1 



Afci = 



rlr2 



(1,., 






V^2 


■ri2- 





12 



} 



(3.12) 



The terms in the first hne of Eqs. p.llj) or ()3.12j) contribute at the 2PN order in the 
conformal factor and they are in agreement with Eq. ()3.5j) . The other terms, proportional 
to «! or 0:2, will not contribute before the 3PN order [because ai and 02 carry a factor 
1/c^ in front]. Therefore, only a small part of the linearized approximation, the one given 
by the first term in 1)3.111) . is necessary in the proof of Theorem ^ As a matter of fact, 
the non-linear perturbation we consider contains much more information than a mere post- 
Newtonian (2PN) expansion. Theorem Q does not constitute a very stringent requirement 
on the non-linear solution of the perturbation equations. 

A point we make by writing the linear solution (T(i) into the primary form ()3.8)1 . involving 
the intermediate functions Hi and Ki, is that the latter functions do enter in the post- 
Newtonian metric at the 3PN order — with the same numerical coefficients as predicted by 
Eq. fl3.8|) . This can be inferred from the expression of the 3PN spatial metric given by Eq. 
(Ill) in Ref. Q|, which contains the particular non-linear potential called X, together with 
the way that the potential X contains the functions Hi and Ki as calculated by Eq. (6.11) 
in jl^ . Thus, although our solution agrees with post-Newtonian calculations up to the 2PN 
order only, it does contain some correct 3PN features. 

C. Quadratic and higher-order approximations 

At the level of the second-order perturbation (oc /3^) the equation reads 

Ao-(2) = ^cr{i)dijSij + djSijdia^i) + Sijdija{i) 

o 



At the next level, cubic-order, the equation will be made of the same terms as in Eq. (|3.13p 
but with the replacements cr(i) a (2) and (T(i), together with many other terms that 
are purely cubic in Sij. And so on for the higher-order equations. 

We shall see in Section IIVI that in order to prove the agreement with the 2PN binary's 
energy (in contrast with the 2PN metric), we need the full information content about the 




(3.13) 
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linearized solution given by Eqs. ()3.1H) - ()3.12p . and also a crucial piece coming from the 
seconc?-orc?er perturbation, solution of Eq. ()3.13p . In the case of non-linear perturbations 
{n > 2), it is in general impossible to find a solution in analytic closed form. Fortunately, 
what we shall need is only to control the expansion of (7(2) at spatial infinity (i.e. in the far 
zone, r +00). And this can be achieved, without disposing of the exact expression of 
(7(2), from the knowledge of the far- zone expansion of the corresponding source-term E(2). 
More generally, the far-zone expansion of the solution o"(„) for any n can be obtained on 
condition that the far-zone expansion of its source has been determined beforehand 
(say, by induction on n). 



The method is issued from the investigation, in Refs. |27l.l28l|. of the multipole expansion 
of the solution of a Poisson equation with non-compact-support source (i.e. whose support 
is M.^). In the present context, the multipole expansion is completely equivalent to the far- 
zone expansion, when r +00. Following Eq. (C.9) in Ref. 3|, we obtain the (formal) 
multipole expansion of (T(„) in the form 



M{a^n)) = Fp|a-i [^(S(„,))] -l-Y^^Llo^ (^1^ ^^rfVl/LS(„)(y)| , (3.14) 

where the calligraphic letter Ai refers to the multipole or equivalently the far- zone expansion. 
The first term in the right side of ()3.14j) represents the effect of integrating "term by term" 
the multipole expansion of the source (this term exists only in the case of non-compact- 
support sources, because the multipole expansion of a compact-support function is zero). 
The second term in Eq. ()3.14j) is parametrized by the "multipole moments" associated with 
the solution, which are given by some definite integrals extending over the non-compact 
support M.^ of S(n)- The symbol FP acts on both terms of (|3^14p . and stands for a certain 
operation of taking the Finite Part, defined in Refs. 13,131 by means of a process of 
complex analytic continuation (with parameter B E C). The finite part has proved to play 
a crucial role, in the case of non-compact-support sources like in order to ensure the 
well-definiteness of the integrals gi ving the multipole moments in ()3.14|) . Notice that Eq. 



fl3.14|) has been proved, in Refs. |27|,|28|, to hold in the case of a regular source [i.e. C°°(M^)]. 



Technical notations in Eq. H3.14(l are for the standard Poisson integral; L = ii • • • if for a multi-index 
composed of £ indices; Ol = • • • c^i^ for the product of £ partial derivatives; ul = Uh ■ ■ ■Uu fo'^ the 
product of £ spatial vectors. We do not write the £ summation symbols over the £ indices composing L. 
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Nevertheless, in the presence of the singular points yi and y2, the formula ()3.14p can still 
be applied, but only in the case of that solution cr(„) for which Act(„) = is satisfied in 
the sense of distributions. 

IV. TOTAL MASS AND INDIVIDUAL BLACK-HOLE MASSES 
A. Asymptotic structure of the solution 

A central result of the present paper concerns the "asymptotics" of our solution — both 
at spatial infinity and in the vicinity of the two particles. 

Theorem 2 (i) The metric 'jij is asymptotically flat at infinity (when r = |x| ^ +00), i.e. 



(a) The two singular points yi and y2 are the images, via some appropriate inversions of 
the radial coordinates: pi ~ 1/ri and p2 ~ l/'"2j of two "asymptotically finite" regions (when 
Pi +00 and P2 — > +00), in the sense that the metric, in coordinates (pi,ni), behaves like 



where Hij denotes a certain function of the angles. 

Theorem |21 says that the solution is composed of an asymptotically fiat universe con- 
nected by continuity, via some Einstein-Rosen-like bridges, to two other regions which are 
asymptotically finite in the sense of Eq. ()4.2|) . Note that, though the metric ()4.2|) . unlike 
the corresponding Brill- Lindquist metric, is not asymptotically fiat in the vicinity of the two 
particles (stricto-sensu), the violation of asymptotic flatness is not very severe, because the 
metric tends toward a constant with respect to pi, and does not involve any divergency "at 
infinity": it remains asymptotically finite — hence the name. On the other hand, what is 
very important is that the real universe, as depicted by this solution, is asymptotically fiat 
at spatial infinity in the usual sense of Eq. 1)4.11) . The global structure described by Theo- 
rem |21 represents an attractive feature, we argue, for considering the solution as a physically 
well-motivated initial state of binary black holes. 




(4.1) 




(4.2) 
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The asymptotic flatness when r ^ +00 follows from the far- zone expansion of the gener- 
ating function ()3.1b|) . which is easily checked to start at Sij = (9 (^) . Using this fact we find 
that the source term for the linearized perturbation — i.e. the right side of Eq. ()H.7|1 — 
behaves like S(i) = O {^)- With the help of Eq. (j3.14p we readily obtain cr(i) = O [^), as 
can also be checked directly with our exact results ()3.11|1 - ()3.12|1 . Then, similarly, we deduce, 
by induction on the non-linear order n [using Eq. ()3.14|) at each step], that cr(„) = O (^) for 
any n. So a = O (^), and the result follows. 

The main point about Theorem |21 is the behaviour of the solution in the vicinity of the 
two particles. When ri — > 0, we find that Sij admits a bounded expansion [i.e. Sij does not 
blow up when ri ^ 0] , of the type 

(3sij = Eijini) + 0{n) , (4.3) 

where Sij is a function of rii = , the direction of approach to the singularity [sij is given 
by Eq. ()A.3a|) ]. From Eq. ()4.3p . the linearized source-term when ri — is like some 
= O (^^y Now by integrating "term by term" the expansion of we get a similar 
expansion, but which starts at the order O . Clearly, the solution o"(i), when ri 0, will 
be composed of the latter expansion, which represents a particular solution, and augmented 
by a possible homogeneous solution, solving the (source-free) Laplace equation. However, 
because our original Poisson equation Ao"(i) = is satisfied in the sense of distributions, 
the only possible homogeneous solution we can add is one of the type "regular at the origin" , 
taking the form of a sum of STF products xl = STF{xi-^ ■ ■ ■ XjJ. Because it is regular when 
ri 0, this homogeneous solution does not modify the leading singular behaviour of cr(i), 
that we have therefore proved to be given by o"(i) = O j . And, again, this can be checked 
with Eqs. (|3.1H) - (j3.12j) . The argument is easily generalized to any non-linear order n (by 
induction on n), thus we conclude that the non-linear perturbation a diverges when ri 0, 
but not faster than a = 0{^). As a result, the expansion of the conformal factor = tp + a 
is of the type 
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where C depends on the unit direction rii but not on ri. From this we deduce that ■jij = \l/^7ij 
behaves dominantly when ri hke 



+£,,(ni) + C(ri) . (4.5) 



ri 

See Eq. flA.4a|) for the expression of C(iii) hnearized order. 

Let us now perform some inversion of the radial coordinate ri. We consider the coordinate 
change, vahd in a neighbourhood of the particle 1, that is defined by x — > (pi, rii), where 

P. ^ ^ . (4.6) 

n 

First we have dx^ = dr\, where = (ri,ni), because the particle is at rest: v\ = 0. Then, 
the coordinate transformation rl — ^ p\ = (pi,ni) involves simply the change of the radial 
variable given by Eq. ()4.6|) . with the angular part rii being unchanged We compute 



dr{ = (^-^^ {6'^ - 1n\n{ + 1n\x^) dp{ , (4.7a) 
^ [6^^ - nin^ ^ . (4.7b) 

(Notice that n\x^ = 0.) We then find that, in the new coordinate system p\ = (pi, rii), the 
metric, say Tij, admits when pi — > +00 the bounded expansion announced in Eq. ()4.2|) . in 
which the quantity iTij (ni) reads 



TT,, = -4nS + 4x'x' + (f^'" - Mn\ + 2n\x') {6'^ - 2n\n{ + 2n\x')eki ■ (4.8) 

It is such that n\n\eij = v^n^nij. [See also Eqs. ()A.5b|) and ()A.6|) .] This completes the 
proof of Theorem El 
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More general coordinate transformations, involving a change of the unit direction ni, are also possible. 
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B. The ADM mass 

Because the metric is asymptotically flat at spatial inflnity (Theorem |2I), the binary's 
total ADM mass is given by the usual surface integral on a topological 2-sphere at inflnity: 

where dS^ is the outward surface element on the surface at inflnity {dS^ = dQr'^rf in the 
case of a coordinate sphere). To be more explicit about M, we recall from Section ll V Al that 
some functions Aij and X of the unit direction n = x/r exist so that, when r —>■ +00, 



p,,. = ^ + O (i) , (4.10a) 

Here Aij is simply linear in while X is in the form of a full non-linear series in 
X = + /3^X(2) + ■ ■ ■ (see also the Appendix). In terms of these functions M is given 

by 



cM , . rdn 



2X + -n'n^Aii 
2 ^ 



(4.11) 



The flrst term is identical to the result in the Brill-Lindquist solution. 

The computation of Eq. ()4.1H1 is quite straightforward at the level of the linearized 
approximation, thanks to the closed-form expression ()3.11|) . However, our aim will be the 
computation of the 2PN energy in Section and we can ascertain beforehand (by counting 
the required powers of G), that among all the terms contributing to the 2PN energy there 
must be one coming from the second-order perturbation, solution of Eq. ()3.13|) . We flnd 
that this particular non-linear term is to be computed only in M (not in the individual 
masses mi, 1712). It comes from that part of o"(2) — or, more precisely, of its leading order 
coefficient X(2) (when r +00) — which does not involve the constants ai and 02, i-e. the 
part which would be the analogue of the flrst term in cr(i) as given by Eq. 1)3.111) . The other 
parts, proportional to ai or 0:2, appear at the 3PN order at least. 

We succeeded in obtaining this non-linear term thanks to the method described by Eq. 
(13.14)1 . First, one can check that the source-term behaves like S(2) = O (pr), so there is 
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no direct contribution, computed "term-by-term" from the expansion of the source-term, to 
the multipole expansion at the order O {1/r): i.e. the first term in the right side of ()3.14j) is 
at least (9(l/r^). The looked-for non-hnear contribution is therefore given directly by the 
value of the "mass monopole" of o"(2), i.e. the integral appearing in the second term of Eq. 
I^rm for £ = 0. So, 



^(2) 



-ji^//Vs,.(y) + o(^) . (4.12) 



Attt 

The integral ()4.12j) can be computed analytically^^, and yields the following contribution to 
the ADM mass: 

where the remainder O (a) indicates that the terms proportional to ai or 0^2 (in this second- 
order perturbation oc are not to be considered for the present calculation. We thereby 
obtain M with sufficient accuracy for controlling the energy at the 2PN order The result 
reads 

The cubic and higher-order perturbations O (/5^) are neglected. See the Appendix for the 
explicit expansion coefficients needed in this computation. 



C. The black- hole masses 



The situation as concerns the two individual masses mi and m2 is less easy than with the 
ADM mass because we dispose only of the notion of "asymptotic finiteness" in the vicinity 
of the particles, described by the fall-off property ()4.2j) . Nevertheless, we wish to find an 
appropriate concept for the black-hole masses. What we shall do is to define nii by the same 

This integral is convergent, thus it is unnecessary to include a finite-part operation FP (for notational 
simplicity we skip the multipole-expansion symbol M). 

This means, by the way, that the relative accuracy on M itself is actually 3PN — because of the rest-mass 
contribution. 
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formula as for the ADM mass M , but using the coordinate system p\ = (pi, rii) in the hmit 
where pi +00. Accordingly we pose 

mi = lim ! dS\(^- ^\ and 1^2, (4.15) 

where Vij{pi, ni) denotes the metric in the coordinates p\, with pi being shown in Eq. ()4.6p . 
and where the outward surface element reads dS\ = dflipln{ in the case of a coordinate 
sphere. 

Because Tij is not asymptotically fiat in the usual sense, due to the term with vTjj in Eq. 
()4.2|1 . the mass defined by the previous integral will be typically unbounded (i.e. mi will 
in general tend toward infinity like pi). Therefore, the definition ()4.15|) does not a priori 
make sense. A possibility would be to discard the divergent part (oc pi) of the mass, and 
thereby to consider only the finite part of the integral in Eq. ()4.15p . i.e. the coefficient of 
the zeroth power of pi in the expansion at infinity. This could represent an appropriate 
postulate for the mass in a general situation [the finite part prescription would be similar 
to the FP present in Eq. ()H.14|l ]. However such a "finite part" process imposed into the 
definition of the mass appears to represent (untill more convincing justification is proposed) 
a somewhat artificial and ad-hoc recipe Gladly enough, we shall not need to invoke any 
ad-hoc finite part because we shall prove that if we restrict ourselves to the computation of 
the binary's energy at the 2PN order, which represents anyway the maximal order at which 
our solution is physically relevant, the masses mi and m2 needed in this computation are 
perfectly well-defined. 

To compute mi and m2, we must control the expansion of the metric when ri — to 
one order beyond Eqs. ()4.3|) - ()4.4|) . Let us pose 



Psij = eij{ni) + ripij(ni) + 0{rl) , (4.16a) 

* = + r7(ni) + C (ri) . (4.16b) 

ri 



The ADM mass, given by the surface term that arises in the Hamihonian, has been obtained even for 
space-times that are not asymptoticahy flat However, it seems difficult to apply this result in the 

present case, notably because of our lack of knowledge of the lapse function N in the vicinity of the 
particles (after radial inversion ri pi of the coordinates). 
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[The expressions of the coefficients are relagated to Eqs. ()A.3|) - ()A.4|) in the Appendix.] It is 
straightforward to derive from this the expression of the metric Tij in the coordinate system 
(pi,ni), taking into account the term ~ l/pi in Eq. (|4.2|) : 



(4.17) 



where Kij = (5''* — 2nfn*^ + 2niX*)(5'-' — 2n[n{ + 2n[x^)nki- Then, by inserting it into Eq. 
fl4.15|) . we obtain 



c2 f dQi 

mi = — hm / —— 

G pi^+oo / 47r 



(4.18) 



Here, njj denotes the trace of the quantity ()4.8jl . and we have used the facts that n\n\'K, 



n\n\eij and together with the useful cancellation of the angular integral 



4:71 



n\n{eij{ni) = , 



(4.19) 



which can be checked from Eq. ()A.4a|) in the Appendix. 

We observe that Eq. ()4.18|) contains a term proportional to pi, with TTjj clS db factor, 
which can and does make the mass to become infinite in the limit pi +oo [indeed, see 
Eq. ()A.7j) ]. As suggested above, such an infinite term could be removed by some procedure 
of taking the finite part (whose meaning would a priori be quite unclear), but we find that 
this infinite term is in fact "negligible" for the present purpose, because Hjj starts only at 
non- linear order in (3, 



- 4x'x' + 4Kx' + nWiX'x')eM = O [p^] . (4.20) 

[See also the expression ()A.6jl .] Therefore we get, at the linearized level, a finite expression 
for the mass (formally): 



rrti 



c" / dVti 



G 



47r 



2r7C(l + n\n{eij) + -("^if^iPij 



(4.21) 



This expression is sufficient to control the 2PN energy. The non-linear corrections O {P^) 
are infinite but will not be needed in this paper. In a sense, they do not belong to the 
"realm" of the present solution, which is limited to the physics at the 2PN approximation: i.e. 
agreement with the 2PN metric in the near- zone, and internal consistency of the asymptotics 
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up to the 2PN order as regards the energy content of the solution (see Section |V|) As for 
the hnearized terms in Eq. ()4.21|) . they are dealt with thanks to the explicit result ()3.1ip . 
and we get 

By setting /3 = into Eq. ()4.22|) . we reproduce the result valid in the case of the Brill- 
Lindquist solution. 



mi 



2c^ai 
G 



V. THE BINARY'S GEOMETROSTATIC ENERGY 

With the masses being defined and computed "geometrically" in Section HVl we get the 
opportunity of an interesting consistency check of our solution, concerning the "geometro- 
static" energy that is associated with those masses. Indeed, in Theorem Q we recovered the 
2PN metric in the near zone, therefore we know the coordinate transformation 5x* = 
between our presently used coordinate system and the harmonic one. The gauge transfor- 
mation vector has been obtained in Eq. ()2.8|) . Therefore it is possible to control, without 
ambiguity — because we determined the coordinate transformation — , the binary's inter- 
acting energy E up to the same 2PN order. If our solution has anything cogent (physically 
speaking), the latter energy should be in complete agreement with the post-Newtonian result 



at the 2PN order known from Refs 



. 0,0, HQ. 



Theorem 3 The binary's interacting energy, deduced from those masses ( |^.^^ and { 4-151 
computed by surface integrals at infinity, i.e. 

E 

— = M-mi-m2, (5.1) 

gives back the energy calculated by standard post- Newtonian methods at the 2PN order (after 
invoking the same coordinate transformation as in Theorem^. 



Recall that our proposal for binary black-hole initial data is to adopt the conformal metric H2.ll) and to 
deduce the conformal factor ^' from it by numerica/ methods. The present analytic investigation [including 
the theoretical definition of the mass 14.1511 ] is aimed at verifying the physical soundness of this proposal. 
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The ADM mass has been previously obtained in the form ()4.14p . while the masses mi 
and 1712 follow from the result ()4.22|) . The corresponding energy E reads then 



E 



c 
G 



4aia2 
ri2 



25(3 \ /3 



+ 



12 



128 



+ 0(a) 



+ O {(3^) 



(5.2) 



As we said above, this formula is accurate enough to get full control of the 2PN approxima- 
tion. We also remarked that at this order the energy is well-defined (i.e. finite). 

Our task is to expand that energy when c +oo, using the facts that ai, a2 = O (^) and 
P = O (^) [actually, we already applied the post-Newtonian approximation when arguing 
that the remainder O {a) is negligible]. Also, we know that (3 is given by Eq. (jH.laj) . We first 
invert Eq. ()4.22|1 . in an iterative post-Newtonian way, so as to obtain ai, a2 with relative 
2PN accuracy. The result is 



Grrii 
~2d^ 



Gm2 G'^m2{5mi + 3m2) 

I ^ _ O 4 I ^ 



2ri2c2 



and 1 



(5.3) 



At the IPN order, we are consistent with the Brill-Lindquist prediction and with what is 
given by Eq. ()3.4|1 . But, quite naturally, we find that Eq. ()5.3|1 differs from the Brill- 
Lindquist result at the 2PN order 

Substituting ()5.3|1 back into ()5.2|1 . we then arrive, after suitable post-Newtonian expan- 
sion, at the expression of the 2PN energy in terms of the two "physical" individual masses: 



Gmim2 G'^mim2{mi + 1712) (j^mim2(mf + 19mim2 + m|) ^ ^ 

It is composed of the standard Newtonian potential energy (for static black holes), aug- 
mented by IPN and 2PN corrections. The result, however, is not yet the one we want to 
prove, because the particles' separation ri2 corresponds to our particular coordinate system. 
If we want to compare E with the post-Newtonian prediction, we must take into account 
the coordinate transformation (in the "near zone" ) that was determined in Theorem ^ 
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In the case of the Brih-Lindquist solution we have 



,BL 



Gmi 



_ Gm2 C^m^jrm +7712) , ( }_ 
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The required link between and the particles' separation Ru in harmonic coordinates 
is readily obtained with the help of Eq. ()2.8|) . which permits to compute the shift of the 
world-lines that is induced by the coordinate transformation We get 

._«_£!M_iI^,Om. (5.5, 



The result we find after replacing Eq. (jS.Sj) into Eq. ()5.4p . and effecting the post-Newtonian 
re-expansion [the replacement is to be made only into the Newtonian term of ()5.4j) ]. is 

Gmim2 G'^mim2{mi + 1712) G^mim2{2m\ + 19mim2 + 2mV\ f ^\ \ 

'' = -^;r^ — i^P +^bJ • ^'-'^ 

Most satisfactorily, we discover that ()5.6p is in complete agreement with the prediction for 
the 2PN energy in harmonic coordinates, calculated in Refs. The expression 

we end up with is the same as given by Eq. (B6) in Ref. js^l — in which of course one must 
set the two particles' velocities to zero. 

The latter agreement is interesting, but we said that it is quite mandatory if our solution 
makes sense. Technically speaking, it necessitated the control of the metric up to second- 
order perturbation theory on the Brill-Lindquist background [cf. the crucial contribution of 
Eq. ()4.1H|1 to the ADM mass]. Furthermore, the result ()5.fi|) can be said to check the global 
character of the solution (notably the asymptotics therein) — because the physical masses 
have been computed by surface integrals at infinity. Alternatively, it shows the relevance at 
2PN order of the definition ()4.15p we adopted for the black-hole individual masses. 

Finally, let us comment that Theorem El tells us something about the physical tenets at 
the basis of the usual post-Newtonian approximation when it is applied to the description 
of black holes. Indeed, the "post-Newtonian" masses mi and 1712, which parametrize the 
post-Newtonian iteration, are introduced as being the coefficients of Dirac delta-functions 
in the Newtonian density of point-like particles Now, we have just seen that in fact 



We discard an infinite "self-term" when considering the coordinate-transformation vector ^i(x) at the 
singular location of the two particles. Thus, from Eq. H2.8|l . which is valid a priori only outside the 
singularities, we obtain the shift vector 

In other words we consider only the finite part of ^i(x) when x — > yi (in Hadamard's sense). 
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these post-Newtonian masses agree (within the present approach at least, i.e. in the time- 
symmetric situation, and up to the 2PN order) with the "geometrostatic" masses which are 
associated with some Einstein-Rosen-hke bridges. We view this as a confirmation that the 
post-Newtonian calculations, which treat formally the compact objects by means of delta- 
function singularities are appropriate to the description of systems of black holes (as long 
as the orbital motion of the black holes can be considered to be "slow" in the post-Newtonian 
sense, i.e. in the so-called inspiralling phase of black-hole binaries). 



Acknowledgment s 

It is a pleasure to thank Gilles Esposito-Farese and Eric Gourgoulhon for discussions and 
remarks. This work was partially supported by the International Academic Grant of the 
Hirosaki University. 



APPENDIX: RELEVANT EXPANSION FORMULAS 

In this Appendix we provide the explicit expressions of the various coefficients in the 
expansion formulas introduced in Section HV1 for the computation of the masses. 

For the ADM mass we need the dominant term, when r +00, in the generating 
function Sij and the solution of the linear-order perturbation for the conformal factor [see 
Eqs. ()4.10|) ]. They are given by 



A 



12 "-12 5 



2ri2 
ri2 I 32 



(n.nia)^ 



+ 



1 «! + ^2 

48 ri2 



(A. la) 
(A.lb) 



where n = x/r and ni2 = (yi — y2)/'''i2 (with n.ni2 denoting the ordinary scalar product). 
The second-order perturbation coefficient was already computed in Eq. 1)4.13^ : 



/9'X(2) 



^12 



256 



+ 0{a) 



(A.2) 
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For instance, the calculation of the 3PN equations of motion of compact binaries reported in Ref. 



H. 
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For the masses mi and m2 we require the expansions when ri of Sij and the conformal 
factor to one order beyond the dominant term [see Eqs. ()4.16|) ]. For Sij we have 





4 


12^12 






^12 '^12 



1 



2 

-ni.ni2 



3 3 
^ni.ni2 + - 



For the conformal factor: 



(A.3a) 
. (A.3b) 



C = «i 



1 + ^1 -^(ni.nis)^ + — ni.ni2 + - 



12 



a2 (3 



1 + — + 



^^12 n2 



24 

1 1 

— Ill .ni2 H 

24 48 



. Oil (I 



o 5 , \ 2 ^ ^ 

' H ni.ni2 ni.ni2 

32^ ' 24 96 



a2 f I 25 

-777111. ni2 - — 

ri2 V 24 48 



(A.4a) 



fO(/3') . (A.4b) 
The quantities and vTjj defined by Eqs. (j4.7bj) and ()4.8|1 read (at hnear order in [3) 



. /5 , 
X = — ni.ni2 - — 
r^2 V 24 

3 



n\2 + (111.1112)^1 



2 I '^f2'^i2 - ^^^^^12 - ^(ni.ni2)nf 
12 \ o o 



(A.5a) 
(A.5b) 



Note that at hnear order TTjj is trace-free. Its trace arises only at second-order in /3, 



/5V 1 

--ni.ni2 + 



12 



6 



144 



:ni.ni2)' - 1 



(A.6) 



At this order the angular average of Tiij is non-zero: 

dQi 5 (3^ 



An 



216 rfa ^ ^ 



(A.7) 
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